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Abstract 

The motion of a particle in 6-dimensional phase space in 
the presence of linear coupling can be written as the sum 
of 3 normal modes. A cubic equation is found for the tune 
of the normal modes, which allows the normal mode tunes 
to be computed from the 6x6 one turn transfer matrix. 
This result is similar to the quadratic equation found for 
the normal mode tunes for the motion of a particle in 4- 
dimensional phase space. These results are useful in track- 
ing programs where the one turn transfer matrix can be 
^f^computed by multiplying the transfer matrices of each el- 



0\ 



ement of the lattice. The tunes of the 3 normal modes for 
motion in 6-dimensional phase space can then be found by 
^solving the cubic equation. Explicit solutions of the cubic 
^equation for the tune are given in terms of the elements of 
K-jthe 6x6 one turn transfer matrix. 

I. INTRODUCTION 

The motion of a particle in 6-dimensional phase space in 
^the presence of linear coupling can be written as the sum 
V^Qof 3 normal mode. A cubic equation is found for the tune 
f^of the normal modes, which allows the normal mode tune 
lOto be computed from the 6x6 one turn transfer matrix. 
This result is similar to the quadratic equation [1] found 
CZ2for the normal mode tune for the motion of a particle in 4- 
^^limensional phase space. These results are useful in track- 
£^ing programs where the one turn transfer matrix can be 
I computed by multiplying the transfer matrices of each el- 
^ement of the lattice. The tune of the 3 normal modes for 
^motion in 6-dimensional phase space can then be found by 
solving the cubic equation. Explicit solutions of the cubic 
equation for the tune are given in terms of the elements of 
the 6x6 one turn transfer matrix. 

II. A CUBIC EQUATION FOR NORMAL 
MODE TUNES 

The particle coordinates are assumed to be x, p x , y, p y , 
z, p z . Then the linear motion of the particle about some 
central orbit can be described by a 6 x 6 transfer matrix 
T(s, s )- 



x(s) 



T(s, so)x(so) 

X 
Px 

y 

Py 

Z 

Pz 
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(i) 



Note that the symbol x is used to denote the column vec- 
tor x and also the first element of this column vector. The 
meaning of x should be clear from the context. This nota- 
tion is used in several places in this paper. 

It is assumed that the equations of motions can be de- 
rived from a hamiltonian, and the matrix T is symplectic: 



TT 
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ST S 



(2) 



T is the transpose of T. 

I is the 6x6 identity matrix, and S can be written in 
terms of 2 x 2 matrices as 



(3) 
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The 2x2 matrix, also called S, is given by 

s _, 1 

" " 1 i) 



(4) 
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It is convenient to also write T in terms of 2 x 2 matrices 



(5) 



One may note that the 2x2 matrices of T are (T)ij = 
Tji. The eigenfunctions of the one turn transfer matrix 
T(s + L,s), where L is the period of the forces acting on 
the particle, obey 

T(s + L, s)x = Xx 

Because T is symplectic then if A is an eigenvalue, 1/A 
is also an eigenvalue [1,2] and the 6 eigenvalue can be ar- 
ranged in 3 pairs of A and 1/A. For stable motion, |A| = 1 
and A can be written as 



A = exp(i / u) 
jj, = 2ttv, 



(6) 



where it is assumed that the period is one turn. 

To find the tunes of the normal modes, one has to find the 
eigenvalues of T, A, which are given by |T — A7| = 0, where 
|T| indicates the determinant of T. It is more convenient 
to find the eigenvalues of the matrix C defined by 



C=i(T + T) 



(7) 



The eigenvalues of C, A, are 



A 
A 



1 
2 

COS jJL 



(A + l/A), 



and the A are determined by 



|C- A7| = 



(8) 



(9) 



It is convenient to write C in terms of 2 x 2 matrices. 
Note that the 2x2 elements of C are 



Cij — (Tij + Tji)/2 



(10) 



and thus 



1 



Tr(Tii)/ 



We can then write C as 



C 



hi C12 
C21 t-jl 
C31 C32 



C13 
C23 
t 3 I 



(11) 



(12) 



1 



Tr (T z ' z ' ) 



To evaluate |C — A7|, C— A7 will be multiplied from the 
left by a sequence of matrices, each with determinant 1, 
to produce a matrix which is upper triangular; that is, 
the elements of this matrix below the diagonal are all zero. 
This triangular matrix has the same determinant as C— A7, 
and its determinant can be found by multiplying all the 
diagonal elements of this matrix. 

First, multiply C— A7 from the left by the matrix 



1 _ 

-C 21 /h 1 
1 



(13) 



A 



This matrix is designed to eliminate C21 and the multipli- 
cation gives the result 



ill _ C12 _ C13 
(t 2 -|C 12 |/ti)7 C 23 - C 21 C 13 /ti 

C31 C32 t3 



(14) 



The result C12C12 = | Ci 2 1 has been used. Note that C21 
Ci2- Now multiply the matrix Eq. (14) by the matrix 



1 


-C31/I1 





1 
1 



(15) 



which replaces the third row in Eq. (14) by the third row 
plus — C31/I1 times the first row giving 



ill _ C12 _ C13 
(? 2 -|C 12 |/? 1 )7 C 23 - C 21 C 13 /?i 
C 32 - C 31 C 12 /h (t 3 - \C 13 \/h)I 



(16) 



One more matrix multiplication is required to reduce C— A7 
to an upper triangular matrix. Multiply matrix Eq. (16) 
by the matrix 



1 

1_ _ _ 

-(C 3 2-C 3 lCi2/<l)/(<2-|Ci2|<l) 



(17) 



This gives 



ill _ C12 _ C13 

(h - \c 12 \/h)i c 23 - c 21 c 13 /li 
7 3 7 



(18) 



7 3 = [t 3 - \c 13 \/h - |c 23 - c 21 c 13 /ii|/(i 2 - IC12IA)] 

The matrix Eq. (18) is a triangular matrix with zero el- 
ements below the diagonal and its determinant, and thus 
|C — A7|, can be computed by multiplying the diagonal 
elements which gives 



|C- AI\ 1/2 
|C- A7I 1 / 2 
|C- A7I 1 / 2 



*l(*2-|Cl2|/*l)[*3-|Cl3|/*l 
-|C23-C 2 lCi 3 /*l|/(*2- |Cl 2 |/*l)] 
h [(*2 - |Ci 2 |/*l)(*3 - |Ci 3 |/*l) 

-IC23-C21CWI1I] (19) 

(1A)[(^2-|C 12 |)(M3-|C 13 |) 
— |C23^1 — C21C13I] 



To simplify Eq. (19), note the following result 

|C23?1 — C2lCi3|7 = (C23?l — C2lCi3)(C23?l ~ C13C21) 

= [|C 23 |*? - T r (C23C 13 C 2 i)ti + |C 2 i||Ci 3 |] / (20) 
which gives 

IC-A7I 1 / 2 = M2<3-<2|Ci 3 |-<3|Ci2|-<l|C23|+T r (C23C 3 lCi 

Finally, the normal modes tunes, which are given by |C — 
A7| = are determined by the cubic equation 

(A-t 1 )(A-t 2 )(A-t 3 )-(A-t 1 )|C 23 |-(A-t 2 )|C3i|-(A-t 3 )|C 

— T r (Ci 2 C 2 3C3i) = 
A = cos /i, 11 = 2tti> 

t'2 



1 



tl = T^ll), 



^T r (t 22 ), * 3 = iT r (T 3 3) (22) 



C 



The Tij are the 2x2 elements of T and Cij are the 2x2 
elements of C. 

Equation (22) is a cubic equation, and the 3 roots of this 
equation A; = cos jj,i,i = 1,3 will give the 3 tunes of the 



normal modes, Vi = Hi/2ir. If Hi is real the mode is stable; 
if n has an imaginary part then the mode may be unstable. 

To get a result for 4-dimensional coupled motion, we put 
C 23 = C31 = and Eq. (22) then gives 

(A-*i)(A-* 2 )- IC12I = 

C12 = ^(Ti2+T 21 ), h = ^T r (Tn), h = ^ r (T 22 ) (23) 

Eq. (23) is the known result [1,3] for the tunes of the 2 
normal modes for coupled motion in 4-dimensional phase 
space. 

In the 4-dimensional case, Eq. (23) can be solved to find 
the two cos Hi of the normal modes. This gives the result 
for cos h 

cos// = i(/ 1+ / 2 )±i[(/ 1 -/ 2 ) 2 + 4|C 12 |] 1/2 

h = ^T r (T n ), *2 = |t,.(T 22 ) (24) 
IC12I = |Ti 2 + T 2 i|/4 

where T 8 j are the 2x2 elements of the one turn transfer 
matrix. If one computes T by multiplying the transfer 
matrices of all the elements in the accelerator, then one 
can find cos /«, and the v values, using Eq. (24). 

Equation (24) is useful in tracking programs for finding 
the normal mode tunes, from the one turn transfer matrix, 
for coupled motion in 4-dimensional phase space. 

A similar result can be found for coupled motion in 6- 
dimensional phase space by solving the cubic equation, Eq. 
(22). 

Equation (22) can be written as 

A 3 + 

a 2 = 

ai = 

a = 

ti = 
dj = 

The solutions of Eq. (25) when the 3 roots are all real, can 
be written as [4] 

cos m = (t 1 +t 2 +t 3 )/3 + 2(r 2 + b 2 ) 1 ' 6 cos(a/3 + S i ) 

Si = 0,27r/3,-27r/3 

b = |r 2 + «/ 3 | 1/2 

r = -(aia 2 - 3a ) - — a 3 (26) 

1 1 2 
q = — ai a., 

H 3 9 2 

tana = b/r 



The condition for the 3 roots to be real is [4] 

q 3 + r 2 < (27) 

If (q 3 + r 2 ) > 0, then 2 roots are imaginary and the motion 
may be unstable. In order to get q 3 + r 2 < and stable 
motion, one has to have q < 0. 

Equation (26) can be used in a tracking program to find 
the 3 normal mode tunes from the one turn transfer ma- 
trix. The one turn transfer matrix can be computed by 
multiplying the transfer matrices for each element in the 
lattice. 
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a 2 A + aiA + a = 

-(h+t 2 + t 3 ) 

tlt'2 + t'2t 3 + t 3 tl ~ |Cl2| " |C 23 | " |C 3 l| 
— ^1^3 — T r (Ci 2 C 2 3C3i) +ti I C23 I + ^2|C3i I + ^3 |Ci2 | 

^T r (T n ), h = ^T r (T 22 ), / 3 = iT r (T 3 3) 
1 



2 C^'i + 



(25) 



Note that the 3 values of Si, Si = 0, 27r/3, = 27r/3 will give 
the 3 roots from Eq. (26). 



